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Abstract. We consider the expulsion of proton fluxoids along neutron vortices from the superfluid/superconductive 
core of neutron star with weak (B < 10 10 G) magnetic field. The velocity of fluxoids is calculated from the balance 
of buoyancy, drag and crustal forces. We show, that the proton fluxoids can leave the superfluid core sliding along 
the neutron vortices on a timescale of about 10 7 years. An alternative possibility is that fluxoids are aligned with 
the vortices on the same timescale. As the result, non-aligned surface magnetic fields of millisecond pulsars can 
be sustained for > 10 9 years only in case of a comparable dissipation timescale of the currents in the neutron 
star crust. This defines upper limits of the impurity concentration in the neutron star crust: Q < 0.1 if a stiff 
equation of state determines the density profile. 
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1. Introduction 

There are observational evidences that the magnetic fields 
(MFs) of millisecond pulsars (MSPs) are stationary over 
the age of > 10 9 years. The spin-down ages of many 
MSPs and the cooling ages of their white dwarf compan- 
ions exceed 10 9 years (Kulkarni 1986; Danziger et al. 1993; 
Kulkarni et al. 1991; Bell et al. 1995). MSP statistics in- 
dicates also, that the ages of MSPs must be greater than 
10 9 years (Bhattacharya & Srinivasan 1991 and references 
therein). 

We assume that both protons and neutrons are in a su- 
perfluid state in the whole core (Baym et al. 1969). The 
softer the EOS and the more massive the NS, the larger 
is the probability that protons and neutrons in the cen- 
tral region of the core are in a normal state (Page 1998a), 
which implies a completely other core field evolution (see 
e.g. Goldreich & Reisenegger 1992). On the other hand, 
the millisecond pulsars are very old NSs. Their internal 
temperature is expected to be much less than any non- 
zero superfluid transition temperature (Page 2001, private 
communication), favouring the superfluid state in the en- 
tire core. 

The superfluid phase transition takes place early in NS's 
life (Page 1998b) and the MF penetrates the core as an 
array of proton flux tubes (fluxoids), each of them carry- 
ing a quantum of the magnetic flux $o = hc/2e « 2 • 10 -7 



G cm 2 . Inside the fluxoids the MF is as high as B p ~ 10 15 
G, and it decays outwardly exponentially with the London 
penetration depth A ~ 10~ n cm. The core superfluid ro- 
tates forming a discrete array of neutron vortices. The 
kernels of vortices and fluxoids consist of normal matter, 
their radii correspond to the coherence lengths £ n and £ p , 
respectively. 

Because of a superfluid drag effect the neutron vortices 
are also magnetized (Alpar et al. 1984). The MF inside 
each vortex is B n ~ 10 15 G. When a vortex crosses a 
fluxoid, the interaction of their MF results in a potential 
barrier, E m ~ 10 MeV per intersection (Ding et al. 1993, 
DCC hereafter), corresponding to a pinning force F m ~ 
E m /X ~ 10 6 dyn per intersection. 

The interaction energy due to density perturbations in the 
center of a fluxoid is E p ~ 0.1 — 1 MeV per intersection 
and causes a pinning force F p ~ £ p /£ n , which is in the 
same order of magnitude as F m (Sauls 1989). 
The strong interaction between neutron vortices and 
proton fluxoids led many authors to develop models 
which consider the MF evolution together with the NS 
spin-down. Srinivasan et al. (1990) and Jahan-Miri & 
Bhattacharya (1994) calculated the expulsion of the mag- 
netic flux from the core assuming that fluxoids and vor- 
tices always move with the same velocity. DCC and 
Jahan-Miri (2000) calculated the magnetic evolution in 
more realistic models, in which other forces, acting on the 
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fluxoids (buoyancy, drag, tension) were taken into consid- 
eration. Konenkov & Geppert (2000, 2001) (KG00, KG01 
hereafter) extended these models by taking into account 
the back-reaction of the crust ( "crustal" forces) onto the 
rate of flux expulsion. All these authors reported the oc- 
currence of a long-living (~ 10 10 years) low MF com- 
ponent, which could be responsible for a non-decaying 
"residual" MF of MSPs. 

In the present letter we calculate the flux expulsion 
timescale in the core of MSPs when the fluxoids slide along 
(parallel) to the neutron vortices, i.e., when the interpin- 
ning between vortices and fluxoids can be neglected. This 
possibility was first mentioned by Muslimov & Tsygan 
(1985). We show, that this timescale is about 10 7 years. 

2. Description of the model 

Besides the force exerted by the neutron vortices, there 
exist other forces which act onto the fluxoids in the NS 
core. The buoyancy force, acting per unit length of the 
fluxoid, is given by (Muslimov & Tsygan, 1985): 
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where R c is the radius of the NS core and e r is the unit 
vector in radial direction when spherical coordinates are 
used. The buoyancy force acts always radially outward. 
The drag force, which arises due to the scattering of the 
relativistic electrons on the MF of the fluxoid, is given by 
(Harvey et al. 1986; Jones 1987): 
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where n c and E-p are the number density and the Fermi en- 
ergy of electrons, c is the speed of light, and v p the velocity 
of the fluxoid's motion through the core. The drag force 
acts opposite to the velocity of the fluxoid. This equation 
is valid only if collective effects are ignored (see DCC). 
In MSPs the proton fluxoids are indeed expected to be un- 
able to cut through the neutron vortices, and the motion 
of arbitrarily oriented fluxoids in the direction perpendic- 
ular to the rotational axis occurs only on the spin-down 
timescale of MSPs. However, the displacement of the flux- 
oid parallel to the neutron vortex does not change the pin- 
ning energy, thus the movement of the fluxoids parallel to 
the vortices is not restricted by the pinning force. There 
is a component of the buoyancy force, /b||j which is not 
compensated locally by the pinning force (Fig. 1). Under 
the action of this non-compensated component, proton 
fluxoids can be either expelled from the core (line 1 in 
Fig. 1), or aligned with the neutron vortices (line 2 and 2' 
in Fig. 1). 

The total flux of the Poynting vector through the crust- 
core interface is equal to the power of forces, acting upon 
the fluxoids in the core (for details see KG00): 
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Fig. 1. The superfluid core of a NS rotates with the angu- 
lar velocity f2 by forming array of neutron vortices (shown 
in the figure by thin vertical lines). Fluxoids are shown by 
dashed lines. The fluxoid 1, which does not intersect the 
equatorial plane of the NS, will be expelled from the core 
along the vortices under the action of noncompensated 
component of the buoyancy force /b||j while its compo- 
nent /bj_ will be compensated by the vortex acting force. 
The fluxoid 2, which intersects the equatorial plane, will 
become aligned with the vortices into position 2'. The 
timescales of both processes are the same. 



where the integral on the l.h.s. is taken over the length of 
a fluxoid, the summation runs over all fluxoids, and the 
integration on the r.h.s. is performed over the crust-core 
interface. For simplicity we assume that /b| = /b (this 
simplification leads to an underestimation of the time of 
expulsion in the order of unity) , that the number density 
of fluxoids is uniform, and calculate all quantities at the 
density of the crust-core interface. Thus we estimated the 
integral on the l.h.s. of Eq. (|) as 4/3 • (/ b + f v )v p N p R c 
(see KG00), where iV p = t:B c R 2 /(& is the total number of 
fluxoids and B c the strength of the mean core MF. Note, 
that Eq. (0) contains no force exerted by neutron vortices, 
and describes the motion of fluxoids along vortices. 

We assume the MF to be axisymmetric, poloidal, 
and dipolar outside the NS. Thus, one can introduce 
the vector-potential A = (0,0, A$), where = 
S(r,t)sin(6)/r — B R 2 s(r,t)sm(9)/r, s is Stokes' stream 
function, normalized to BqR 2 , Bq and R are the initial 
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Fig. 2. The top panels show the evolution of MF strengths in the core (B c ), at the surface (B), and at magnetic 
equator at the crust-core interface (Bg). The bottom panels show forces, acting on the fluxoids in the core of MSP. 
Left, middle and right columns correspond to Q = 0.1,0.01 and 0.001. 



surface MF strength and radius of the NS, respectively. 
The integral on the r.h.s. of Eq. (||) is given by 
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We introduce the total core forces Fb jV = fb,v • 4i? c /3 • N p , 
and the "crustal" force, which arises when the roots of the 
fluxoids are moved along the crust-core interface: 



Eq. 
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can be rewritten as 



dr 



F h + F v (v p ) + F cmst (v p ) = 0. 



(5) 
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The "crustal" force i^crust depends on v p in a rather com- 
plicated way, since both ds(R c , t)/dt and ds(R c , t)/dr de- 
pend on v p . The evolution of the homogeneous MF within 
the core is governed by v p too (for details see KG00). F crus t 
might be found by solving the induction equation in the 
NS crust {R c <r < R): 
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where a is the conductivity of NS crust. The crustal con- 
ductivity is determined by collisions of the electrons on 
impurities and phonons. For the phonon conductivity we 
use the numerical data given by Itoh et al. (1993), for the 
impurity conductivity we apply the analytical expression 
derived by Yakovlev & Urpin (1980). The magnitude of 
the impurity conductivity is determined by the impurity 
parameter Q, whose value is highly uncertain and has to 
be varied. The impurity conductivity dominates in the in- 
ner crust, where currents, generated by the flux expulsion, 
are expected to be localized. 



The boundary conditions are given by s(R,t)/R = 
—ds(R, t)/dr at the surface which joins the MF inside the 
NS with the dipolar field in the vacuum outside, and by 
Eq. (|J) at the crust-core interface. The initial condition 
is set by Bo and the initial s-profile in the NS. 

3. Results 

In Fig. 2 we show the evolution of the MF strengths and 
of the forces acting on the fluxoids for Bq = B c q = 10 9 
G and an initial spin period Pq = 0.001 s. The results 
in the left, middle and right column are obtained for an 
impurity parameter Q — 0.1,0.01 and 0.001, respectively. 
The computations were performed for a NS model, based 
on the Friedman-Pandharipande equation of state (EOS), 
with NS mass M = 1.4M Q , radius R = 10.4 km, and 
thickness of the crust AR = 934 m (Van Riper 1991). 
The velocity of the fluxoids remains almost constant (ss 
7 • 10~ 9 cm/s) during the MSP evolution for all values of 
Q considered here. This is because in Eq. (0) i^ crus t can be 
neglected in comparison with and F v during the whole 
evolution. It is seen from Fig. 2, that F crust < O.Ol^b, F v ) 
and the fluxoid velocity is determined only by the balance 
of buoyancy and drag forces. The expulsion timescale of 
the MF from the core, R c /v p , is < 10 7 years for all val- 
ues of Q in the crust, whereas the decay timescale of the 
surface MF is dependent on Q, and can be estimated by 
10 7 /<5 years for the given EOS. Note, that the timescale 
of the alignment must be in the same order of magnitude 
as the expulsion timescale. 

An interesting feature is the generation of a strong 
^-component of the MF at the crust-core boundary. 
Depending on Q it may exceed the surface field strength 
by up to three orders of magnitude when the flux expul- 
sion proceeds. This effect reflects the induction of currents 



4 



D. Konenkov, U. Geppert: Residual Magnetic Fields in Millisecond Pulsars 



close to the crust-core interface, which are the stronger the 
longer the characteristic decay time in the crust is. 
We have performed similar computations for Bq = 10 8 and 
10 10 G and values of the impurity parameter in the interval 
10~ 4 < Q < 0.1. We obtained the same results: the expul- 
sion timescale of the flux from the core is < 10 7 years, 
the timescale of decay of the surface MF is 10 7 /Q years. 
That means that the expulsion timescale does not depend 
on the MF strength and Q if B a < 10 10 G. This is because 
both Ft, and F v are proportional to the core MF strength, 
while i^ crus t is negligible. Our estimate for the flux expul- 
sion timescale in case of MF strengths, characteristic for 
MSPs, is consistent with that of Jones (1987). Note, how- 
ever, that for standard pulsars with B ~ 10 12 — 10 13 G, 
the expulsion timescale is dependent on the MF strength. 
This is because F crust , which is proportional to the -B 2 , is 
then the dominating force, which hampers the expulsion 
of the flux from the core of canonical pulsars (KG00) and 
determines the expulsion timescale, balancing Ft,. Though 
our model simplifies the real magnetic configuration, we 
expect that it gives the correct estimate for the expulsion 
timescale along the rotational axis (10 7 years). In con- 
tradiction, the expulsion across that axis occurs on the 
spin-down timescale, which for MSPs exceeds 10 9 years. 

4. Discussion 

If in MSPs the dissipation timescale of crustal currents 
is less than 10 9 years, their surface MF would vanish. 
Alternatively, the surface MF of MSPs can be aligned 
with their rotational axis on timescales of < 10 9 years. 
Both effects would make MSPs unobservable as radiopul- 
sars. The observed non-aligned long-living surface MFs of 
MSPs can be sustained only if the dissipation timescale of 
currents in the NS crust is > 10 9 years. Because of large 
differences in the thickness of the crust, for a soft EOS 
the decay timescale of the crustal MF is roughly given 
by 10 6 /Q years, while for stiff EOS that relation reads 
10 8 /<3 years (Urpin & Konenkov, 1997). Thus, for a soft 
EOS and a thin crust the impurity parameter must be as 
low as Q < 1CP 3 , while for a stiff EOS with a thick crust 
Q can be much higher (Q < 0.1). 

The models of flux expulsion, developed by DCC, KG00, 
KG01, and Jahan-Miri (2000) for standard pulsars, give an 
inadequate low velocity of fluxoids at the late evolutionary 
stages, resulting in the appearance of long-living residual 
MF components. The main reason is that the equations of 
balance of forces, which determine the velocity of fluxoids 
(see Eq. (31) in DCC, Eq. (11) in Jahan-Miri 2000, Eq. (6) 
in KG00), arc written in scalar form, whereas the forces 
are vectors. These equations describe quantitatively only 
the evolution of the component of the fluxoid velocity di- 
rected perpendicular to the vortices. We now rediscuss the 
papers KG00, KG01 and show, that this is a good approx- 
imation for two evolutionary stages. I) In relatively young 
(t < 10 4 years) standard radiopulsars with B ~ 10 12 — 10 13 
G this velocity component is determined by the balance of 
vortex acting force F D and the drag force F v , and is about 



10~ 8 - 10~ 7 cm/s (KG00). The force exerted by fast out- 
ward moving neutron vortices exceeds the buoyancy force. 
Thus we expect that the component of v p perpendicular 
to the rotational axis would also exceed the parallel com- 
ponent of v p . II) Later, when F n <C (Ft,, F crust ), Ft, is bal- 
anced by -Fcmst thus determining the expulsion timescale. 
Then, the fluxoids are almost not affected by the vortices 
when they are moving towards the crust-core interface 
across them. At this stage the bulk of the flux is expelled 
from the core, and we conclude that the timescale of ex- 
pulsion was estimated correctly. 

When almost all magnetic flux is expelled from the core, 
the vortex acting force is balanced by the perpendicular 
component of buoyancy, thus hampering further flux ex- 
pulsion across the vortices. However, expulsion may pro- 
ceed along the vortices on a timescale of ~ 10 7 years. The 
appearance of the residual MF in the above cited papers 
seems to be a qualitatively wrong result. We do not ex- 
pect any unaligned magnetic flux to remain in the cores of 
old (10 9 - 10 10 years) NSs with low (B < 10 10 G) surface 
MFs. 

5. Conclusion 

The basic idea of this paper is that the motion of proton 
fluxoids parallel to the neutron vortices is not restricted 
by the pinning force, since the displacement of (even 
pinned) fluxoid parallel to the vortex does not change the 
pinning energy. This approach is basically different from 
those which consider the occurrence of residual fields by 
the balance of buoyancy and vortex acting forces at the 
late stage of NS evolution. 

The velocity of fluxoids in the superfluid core of the NS 
with low MF (B < 10 10 G) is calculated. It has been 
shown, that the fluxoids escape from the core or become 
aligned with the neutron vortices on a timescale of 10 7 
years, determined by the balance of the buoyancy and 
drag forces. This timescale is independent of the MF 
strength and of the crustal conductivity. The surface 
MF of the MSPs can be maintained on the timescale of 
> 10 9 years only if the dissipation timescale of the crustal 
currents is also > 10 9 years. Therefore, an estimate of 
the impurity parameter Q is possible. Once the NS core 
matter behaves according a soft EOS this gives a limit 
of Q < 0.001, for a stiff EOS Q < 0.1, counteracting the 
slower field decay in a thicker NS crust. Because of the 
relatively weak MF in MSPs neither the spin-down, nor 
the "crustal" force play an important role for the flux 
expulsion. 

Many uncertainties affect the estimation of Q, both for the 
crust of isolated and accreting neutron stars. However, it 
is very unlikely that the impurity parameter is as small as 
10~ 3 . Dc Blasio (1998, 2000) argues that Q is relatively 
large in that layers of the crust of isolated NSs, where 
the chemical composition changes; in average the impu- 
rity parameter should be larger than 10~ 3 . The calcula- 
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tions of the chemical composition and impurity content 
of the crusts of accreting NSs show that Q can be even 
as high as ~ 100 (Schatz ct al. 1999). Since it is very 
likely that MSPs went through a phase of intense accre- 
tion, this indicates - under reserve that the suppositions 
of our model are fulfilled - that at least not too soft EOS 
should describe the state of NS core matter. Such a conclu- 
sion has been drawn also from the observation of kilohertz 
quasi-periodic oscillations in several low mass X-ray bina- 
ries (Kluzniak 1998) as well as from the consideration of 
the magneto-rotational and thermal evolution of isolated 
NSs with crustal MFs (Urpin & Konenkov 1997; Page et 
al. 2000). 

Clearly, the scenario will be qualitatively changed when 
a portion of the core's volume remain in the normal state 
or the forces acting on the fluxoid (Eq. (2)) have to be 
modificatcd due to collective effects. 
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